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In this paper, virtual reality techniques are 
used to define an intuitive interface to a 
nanoscale manipulation device. This 
device utilizes optical methods to focus 
laser Hght to trap and reposition nano-to- 
microscopic particles. The underlying 
physics are simulated by the use of 
Lagrange mechanics. A unique control 
method for the manipulation of the parti- 
cles is also provided. The user can natural- 
ly grab and steer the particles. Behind the 
scene, a complex computation is per- 
formed to find the new location of the 
potential field induced by the laser beam 



that would move the particles accordingly. 
Haptic feedback is used to constrain the 
steering motion within the physical capa- 
bility of the potential field. 

Key words: nanoscale assembly; nan- 
otechnology; optical tweezers; virtual real- 
ity. 

Accepted: June 3, 2003 

Available online: http://www.nist.gov/jres 



1. Introduction 

The development of nanotechnology [1] has generat- 
ed a great deal of scientific interest, but the promise of 
nanotechnology will be realized only by manufacturing 
commercial products from assemblies of atoms, mole- 
cules and nanoscale components. While nanoscale 
measurement and fabrication techniques are advancing 
rapidly few effective techniques for manipulation of 
nanoscale objects are available. Manipulation of 
nanoscale objects can be achieved using a scanning 
probe microscope, such as an AFM (Atomic Force 
Microscope) or an STM (Scanning Tunneling 
Microscope). Guthold et al. [2] performed rolling and 
sliding of carbon nanotubes using an AFM. Eigler and 
Scheweizer [3] used an STM to position individual 
xenon atoms on a single-crystal nickel surface, but 
these tools do not yet offer very meticulous control of 
objects (e.g., the ability to grasp and release). Laser 
beams can also be used to trap and manipulate small 
particles. A laser apparatus, called OT (optical tweez- 



ers) [4], provides the user with a noncontact method for 
manipulating objects that can be applied to viruses, 
bacteria, living cells and synthetic micro and nanoscale 
particles. 

Each method has strengths and weaknesses. The 
AFM can be used to manipulate a larger set of objects, 
and the relative forces required to manipulate the parti- 
cles can be more easily derived. The disadvantage is 
that it is generally less accurate than the STM, and 
since the AFM can directly contact the sample, consid- 
erations of sharpness and wear of the tool tip must be 
accounted for. The STM, being non-contact, has little 
tool wear, yet manipulation is limited to conducting and 
semiconductor materials and typically will occur under 
vacuum conditions. The strength of OT is that they are 
well suited for manipulating biological objects and 
refractive and metallic particles, yet it is usually used in 
solution. The AFM and STM are generally limited to 
two dimensions with a very limited third dimension 
(sometimes referred to 2.5 dimensions) whilst OT can 
work in three dimensions. It is also interesting to note 
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that AFM and STM can work both in imaging and 
manipulation mode, but OT has primarily been used for 
manipulation, although 3D imaging has been reported 
[5]. All these techniques require complex driving elec- 
tronics under sophisticated computer control. 

The operation of nanoscale manipulators can be done 
through tele-operation or automatic manipulation. In 
the former approach, a human operator acts as a part of 
the control loop and issues commands by judgments 
based on various forms of sensory cues. The sensory 
cues can be audio, visual and haptic. In the latter 
approach, the computer program makes autonomous 
decisions based on an underlying set of constraints or 
rules to perform the predetermined manipulation sce- 
nario. In developing these constraints and rules, one 
must have some understanding of nanoscale dynamics, 
but this is still an emerging field of research. This is 
why tele-operation is gaining favor for certain user 
interfaces. Tele-operation needs to virtually present the 
nanoscale environment to the operator such that the 
operator can be aware of the target objects and the sur- 
roundings. Various techniques in virtual reality can be 
used to enhance this human-in-the-loop control system 
by letting the user feel immersed in the environment. 
Among these techniques, haptic display plays a vital 
role for the following reason. 

One of the most challenging problems in nanoscale 
manipulation is that real-time image sensing of 
nanoscale artifacts is difficult. The imaging sampling 
time of an AFM and STM can take seconds, or possibly 
minutes, depending on the scan rate and the sampling 
range, and cannot be done simultaneously with the 
manipulation. Thus the only real-time feedback occurs 
with the AFM by way of the force imparted on the 
probe when sensing nanoscale objects (contact or non- 
contact modes) or by the STM current fluctuations that 
can be correlated to height of the tip above the surface. 
The user feedback occurs through linking the AFM 
force or STM height measurement directly to a haptic 
device using an appropriate amplification routine. This 
need for real-time exploration is what attracts 
researchers to develop man-machine interfaces for 
nanoscale manipulation that use haptic display tech- 
nologies. Some of the work is described below. 

Mollis et al. [6] demonstrated that atomic-scale land- 
scapes could be explored and felt with the hand in real 
time by interfacing a haptic feedback device with an 
STM. Their haptic device is based on Lorentz magnet- 
ic levitation. Sitti and Hashimoto [7] used an AFM with 
a one degree of freedom haptic device to sense 
nanoscale forces. They demonstrated two-dimensional 
positioning of micrometer sized latex particles. 



Guthold et al. [2] used a force feedback device coupled 
with an AFM tip to haptically sense and manipulate 
carbon nanotubes, DNA and viruses. All of these 
research efforts used a graphical display in addition to 
the haptic display. The topographical data used for the 
graphical display of the nanoscale objects were 
acquired by using the probe in the imaging mode that 
was also used in the manipulation mode. The limitation 
of this method is that while manipulating the specimen, 
the graphical display is static and requires additional 
scans to see the result of the manipulation. 

With OT, typical manipulation can be done through 
a two-dimensional pointing device such as a mouse or 
joystick. Usually the manipulation is done in noncon- 
tact mode, so there is no direct force that can be ren- 
dered on a haptic display. Because of the limited abili- 
ty of traditional optical microscopy to resolve nanome- 
ter-scale structures, OT, unlike AFM or STM, must use 
other means and synthesize virtual sensory inputs to the 
operator. This paper presents the work on developing 
an intuitive interface for manipulating particles through 
OT with the use of VR (Virtual Reality). 



2. Approach 

The Lagrange equation that serves as the foundation- 
al model for the trajectory of a particle exposed to a 
potential field is first introduced. The potential field is 
an abstraction of the potential energy field that is 
induced by the laser beam. Note the potential field 
depends on both the laser and the particle. For example, 
if the geometry of the particle changes, the potential 
field is no longer applicable. Then the equation of 
potential field position that can be used to steer a parti- 
cle is derived. Lastly, an approach in constraining the 
steering motion within the physical capability of the 
optical laser instrument through the use of haptic feed- 
back is presented. 

2.1 Modeling of Optical Tweezers Physics 

The time required to model the physics associated 
with the optical tweezers instrument depends on the 
objective of the work and the accuracies required. 
Maxwell's equations can be used at the fine level, but 
geometric optics calculations can give sufficient expla- 
nations for most cases [8]. The approach described in 
this paper goes one step further and assumes that the 
potential energy can be represented by an explicit spa- 
tial function. For the convenience of implementation 
only two-dimensional space and the motions within this 
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space is discussed, but this involves no loss of general- 
ity for a sphere trapped in a cylindrically symmetric 
beam. 

As an example, assume one has a spherical particle 
of radius r, mass m and with a homogeneous material 
distribution. The form of the trapping potential will 
depend on the size and properties of both the beam and 
the particle, as well as the index of refraction of the sur- 
rounding medium. To focus on the visualization of 
trapped particles, a simple potential well with a 
Gaussian distribution is assumed. Assume that due to 
the laser beam, the approximate potential energy of this 
particle can be given by 



-ln(2) 



-^oe 



(1) 



where Po is the well depth, x andj^ are the two axes per- 
pendicular to the laser path, and a and b are the beam 
waist dimensions. The factor 2 in Eq. (1) normalizes 
the potential so that its value at x = a, j = or at 

x = 0,y = b is equal to — - . The potential field given 

in Eq. (1) can be realized by dithering the laser beam or 
through the use of multiple laser beams. Equation (1) 
can be used to describe a spherical potential when a and 
b are the same or a line potential when one of the beam 
waist dimensions is much greater than the other. 
Assuming there is no change in the rotational speed of 
the particle, the kinetic energy of this particle can be 
given by 



In Cartesian coordinates, Eq. (3) expands to 
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where F^ and Fy are external forces to the system. For 
example, viscous forces can be used as the external 
forces. However, for simplicity, they are not considered 
in the calculations. 

Substituting the Tin Eq. (2) and the P in Eq. (1) to 
Eq. (4) leads to 
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Rearranging the above leads to 



r = ^w(x'+/). 



(2) 



From Eqs. (1) and (2), Lagrangian equations of motion 
[9] allows us to predict the behavior of the particle tra- 
jectory. The behavior is expressed as a system of n sec- 
ond-order differential equations. 
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where ^o'^i'*"'^«-i ^^^ %^%r"^^n-i are the general- 
ized coordinates and the time derivatives of these coor- 
dinates. Qo,Qi,"'>Qn-i are generalized forces. 
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Equation (6) can not be solved analytically, and a 
fourth-order Runge-Kutta method for solving the dif- 
ferential equation is used. With the help of Eq. (6), the 
trajectory of a particle can be simulated. What is more 
interesting is that one can solve for the new location of 
the laser trap (center of the potential field) such that it 
would influence the particle to move in a certain direc- 
tion within a specific time lapse. This can be used to 
steer the particle with greater accuracy. 

2.2 Steering the Particle 



Most of the current work on OT moves the laser trap 
position to the desired destination of the particle and 



277 



Volume 108, Number 4, July-August 2003 

Journal of Research of the National Institute of Standards and Technology 



lets the particle drift to the laser position. For light par- 
ticles or strong lasers, there would be little difference 
between the center of the potential field and the parti- 
cle. But for heavy particles with low power lasers, there 
will be noticeable time delay between the two posi- 
tions. It is much preferable to use lower power lasers 
for many reasons other than being more economical. 
For example, strong lasers can damage the exposed 
particle. It is also true that the attraction force is not 



maximum at the center point. By putting the particle in 
the maximum force region, the steering response time 
can be reduced. When the particle needs to be moved 
quickly the particle is put in these regions to reduce the 
response time. 

With the use of the fourth-order Runge-Kutta 
method, the system of differential equations in Eq. (6) 
can be converted to difference equations. 
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where Ax and Ay are the displacements of the particle 
in the x, y directions. Equation (7) computes the dis- 
placement of a particle due to the potential field in At. 
Inversely, with known At and the velocity of the parti- 
cle, one can compute the location of the particle x, y rel- 
ative to the center of the potential field. In other words, 
the center of the potential field relative to the particle 
location that would move the particle in A^ time with 
the displacement of (Ax, Ay) is -x, -y. To solve for x, y 
in Eq. (7) Newton's method [10] is used. The success of 
Newton's method depends on providing it with a suffi- 
ciently good initial guess. The initial guess used is 
obtained in the parameter domain x and y. The domain 
is first divided into square elements. Now for each ele- 
ment, the linear segments of the curve satisfying the 
first part (Ax) of Eq. (7) is computed. The linear seg- 
ments are determined by the signs of Ax at four comer 
points. This technique is a specialization of the 
Marching Cubes algorithm [11] applied to two dimen- 
sions. For each element with segments satisfying Ax, 
the second segment set satisfying Ay is obtained once 
again. At the end, the possible intersection points 
between the first and the second sets are computed. 
These intersection points are refined through Newton's 
method. If the Newton's method does not converge, the 
element is subdivided into smaller sub elements, and 
the procedure is started all over again. 

2.3 Working Envelope 

The working envelope (WE) is defined as the enve- 
lope of displacement vectors a particle can be moved 
by a potential field. The WE is dynamic in that it 
changes due to the velocity of the particle and the time 
lapse until the next control cycle. WE limits the motion 
of the particle. 

The WE is numerically computed. By using the coor- 
dinate values at equally spaced, regular grid points as 
the displacement vectors of a particle with certain 
elapsed time and initial velocity of the particle, solu- 
tions to Eq. (7) are sought. The grid points where solu- 
tions exist form a lump. By connecting the boundary of 
this lump with a loop composed of line segments, one 
gets the WE. Feldman's algorithm [12] is used for com- 
puting the loop pruning the spurious branches. 
Furthermore, the entire interior set of grid points to the 
loop are checked to determine if they have solutions. 
Each point with no solution is recalculated with small- 
er square elements for computing the initial guess point 
to see if that would lead to solutions. In cases where 
this recursive refinement does not eventually lead to a 



solution, an average of neighboring solutions is used as 
the solution. 

2.4 User Interface 

The operator steers the particle with a cursor, which 
is a small spherical ball that is position controlled by 
the stylus of the haptic device. The operator moves the 
cursor to the particle of interest and presses the button 
attached to the stylus and then the particle is rigidly 
glued to the cursor. The system is now in the steering 
mode. The system then shows the WE of the particle in 
question. The WE depends on the time lapse between 
the user control cycles and the velocity of the particle. 
The haptic stylus shows no resistance when the particle 
is moved inside the WE yet encounters a virtual haptic 
wall when the particle is moved outside the WE. A sim- 
ple linear spring model is used to render the wall. 

The ellipse in Fig. 1 (a) represents the WE of the par- 
ticle at the center. When the commanded displacement 
of the particle is within the boundary of the WE, such 
as the two arrows that are inside the ellipse, the opera- 
tor experiences no reaction force, and the particle is 
moved to that position indirectly by the movement of 
the potential field. In Fig. 1 (b), however, the com- 
manded displacement is outside the boundary of the 
WE. In this case the displacement intersects with the 
WE and the particle is moved up only to the intersec- 
tion point. The excess amount of the displacement com- 
mand (initial displacement command subtracted by the 
modified displacement command) is used to determine 
the relative magnitude of the reaction force. This way 
the operator would sense a virtual wall at the WE. 
Notice this way the particle would still move when the 
issuing command is physically incorrect. This maxi- 
mizes the capability of the optical tweezers. 



3. Results 

The results are divided into three subsections. Firstly, 
understanding of the particle trajectory is gained with 
various experiments. Secondly, the numerical results in 
computing the laser position and WE is given. Lastly, 
the virtual environment is presented. 

3.1 Characteristics of the Particle Movement 

The trajectories of particles in the potential field can 
be computed by Eq. (7). To gain insight into how the 
particles would move in the potential field, the trajecto- 
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Fig. 2. Trajectories of particles in optical potential field (a) one time 
step of 1 fxs (b) ten time steps of 0.1 |lis. 



Fig. 1. Displacement of the particle (a) within WE (Working 
Envelope) (b) outside WE (Working Envelope). 

ries of particles placed in grid points starting from rest 
is plotted. The result is illustrated in Fig. 2 (a). 
Following constants are used: a = 600 nm, b = 300 nm, 



Pq=1,0x 10^ nJ, m = 1.0 X 10"'^ g, x = nm/|Lis, y = 
nni/|Lis and A^ = 1.0 |Lis. The grid spacing in Fig. 2 is 160 
nm. To examine the nonlinear nature of the trajectory, 
the experiment is tried with multiple time steps. The 
result of using 10 time steps totaling one |Lis is given in 
Fig. 2 (b). It can be seen that a more accurate result can 
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be achieved by using multiple time steps, but it requires 
more computation time. As mentioned before, fourth- 
order Runge-Kutta with a fixed time step was used for 
the numerical integration. The calculated results agreed 
with the results from the commercial numerical solver 
(SIMULINK^). A varied time step method was also 
tested in SIMULINK giving agreeing results. 

The next question was how far the particles would 
travel in the given potential field. Figure 3 (a) is a plot 
of the displacement distance as a function of the parti- 
cle location from the center of the potential field. 



Figure 3 (b) is the plot of iso-contours having the same 
displacement distance. The legend of the distances is 
shown at the top right comer. Notice that two spots 
symmetrical to the horizontal axis have maximal dis- 
placement. These spots can be important as they dictate 
the displacement limit of a given laser force field. It is 
also important to note that the center spot of the laser 
beam shows low displacement distance. The locus of 
maximal displacements in radial directions is given in 
Fig. 4 (a). The directions were sampled uniformly in 
increments of 3.6°. 
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Fig. 3. Plot of displacement distances (a) magnitude (b) iso-contours of equal distances. 



Certain commercial equipment, instruments, or materials are iden- 
tified in this paper to foster understanding. Such identification does 
not imply recommendation or endorsement by the National Institute 
of Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the pur- 
pose. 
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In Fig. 3, the particle was initially at rest, the kinetic 
energy being zero. Since the potential energy is higher 
on the periphery than at the center, the combined total 
energy is higher on the periphery. This means that each 
(x, y) location in Fig. 3 depicts a particle in a different 
total energy state. One can obtain a different plot when 
one enforces constant total energy. Fig. 5 (a) is the 
coordinate system used for the plot in Fig. 5. The 
ellipse denotes the constant potential energy curve. The 
small black circle lying on the ellipse represents a par- 
ticle with a velocity at a certain angle ("velocity angle") 
to the horizontal axis. The position on the curve is 
determined by the "position angle." Angles in Fig. 5 (b) 
and (c) are shown in radians. The particle displacement 
distance is plotted by varying the particle position and 
its velocity direction on a constant potential energy 
(1.0x10"^ nJ) curve with constant kinetic energy 
(5.0 X 10"^ nJ). The result is shown in Fig. 5 (b). Figure. 
5 (c) is the plot of iso-contours having the same dis- 
placement distance. 
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Fig. 4. Locus of maximum displacements in radial directions (a) 
particle initially at rest (b) particle initially moving. 



In realistic situations the particles would have veloc- 
ities. The effect of velocities on the particle trajectory is 
studied. In Fig. 4 (b), the following constants are used: 
i=100 nm/|Lis, j=100 nm/|Lis. Notice the locus is 
changed due to the initial velocity of the particle. 



3.2 Computing the Laser Position 

Now comes the most important part, computing the 
center of the potential field such that a particle at a cer- 
tain initial velocity would displace to a required posi- 
tion in a given amount of time. The solution of Eq. (7) 
can be very involved due to the highly nonlinear nature 
of Ax and Ay, especially when At is large. The plots of 
Ax and Ay are given in Fig. 6 (a) and (b), respectively. 
Following parameters are used: A/= 1.3 s, x=100 
nni/|Lis, j = 100 nni/|Lis. For visual clarity the axes are 
moved to the bottom left. The origin is at the center. In 
Fig. 7, the two small loops aligned vertically at the cen- 
ter are the locus of points satisfying Aj^ = -800 nm. 
Similarly, the biggest loop is the locus of points satisfy- 
ing Ax = -160 nm. The two intersection points shown 
in Fig. 7 are the two solution points satisfying both con- 
ditions. If this is extended spatially, and the solutions 
are sought, one can compute discrete points where the 
solutions exist or do not exist. For example, when one 
uses the coordinates of regular grid points as Ax, Aj^ and 
use small dots to represent points where there are no 
solutions and larger dots for points where there are, one 
can get Fig. 8 (a). The initial guess points were 
obtained with the grid distance equally spaced as 
shown in Fig. 8 (a). The computation took 17 s on a 
Pentium PC. Notice that the envelope of solvable grid 
points is connected with line segments. The connected 
segments form the WE. If desired, one can get more 
accurate results by using a finer grid distance for the 
initial guess points. Figure 8 (b) was obtained with four 
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Fig. 6. Solution surface (a) plot of Ax (b) plot of A>'. 



times finer resolution and took 288 s. Furthermore, Fig. 
8 (c) was obtained with sixteen times finer resolution 
and took 4591 s. Notice there is very little difference 
between Fig. 8 (b) and Fig. 8 (c). Formal formulation 
for optimal grid spacing for the initial guess point is not 
easy to do. Therefore the results as a function of differ- 
ent resolutions are examined and the one where the 
increase in the resolution did not significantly improve 
the solution was chosen. 

3.3 The Virtual Environment 

Figure 9 (a) illustrates the initial screen of the virtu- 
al environment (VE). The descriptions of the graphical 
objects are as follows. The solid sphere at the center is 



the particle. It is being illuminated by a dithered laser 
beam shooting in the direction of the plane. This is 
shown as the ellipse in Fig. 9 (a). The small sphere to 
the left of the particle is the cursor for the operator. 
When the haptic stylus is moved, the cursor follows 
accordingly. The particle can be moved only after it has 
been grabbed. To grab the particle, the operator moves 
the cursor to either touch or be inside the particle and 
presses the button attached to the stylus. To release the 
particle the operator simply releases the button. The 
diameter of the particle is one jim and the grid spacing 
distance is 2.9 jim. The physical parameters used for 
the potential field are a = 600 nm, b = 300 nm,pQ = 1.0 
xlO^nJ, m=1.0xlO-''g. 
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Fig. 7. Solution points. 

Figure 9 (b) illustrates the instant snap shot when the 
particle is steered to the right and the displacement 
amount exceeds the physical capability of the potential 
field. The virtual environment reacts by reacting with a 
force that is in the direction counter to the steering 
motion. The amount of force is proportional to the 
exceeding distance from the WE. The PHANTOM hap- 
tic device was used for the haptic output and steering. 
This was attached to a Windows NT with a Pentium 
processor and 256 MB of RAM. The descriptions of the 
graphical objects are as follows. The arrow direction 
and length shows the direction and magnitude of the 
reaction force applied to the haptic stylus. The slender 
rod with elliptical cross section denotes the dithered 
laser beam. Finally the closed loop shows the WE. 
Because the WE is too small, it was drawn with a scale 
factor of 100. The WE and its interior solution points 
are pre-computed before the execution. The bounding 
box region with range of -100 nm to 100 nm on each 
axis was used to compute the WE. The cached data 
depends upon the elapsed time and the initial velocity 
of the particle. For simplicity and the limitation of the 
computer memory, a constant cycle time of 0.2 |Lis was 
used. The initial velocity of the particle at each axis was 
varied from -80 nm/|Lis to 80 nm/|Lis in increments of 
5.33 nm/|Lis. To introduce a viscous effect and also to 
narrow down the range of velocity of the particle to be 
considered, the velocity was reduced to one tenth of the 
preceding cycle at each cycle. Note that a more formal 
way to have done this would be to have modeled the 
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Fig. 8. Working envelope (a) 20x20 grid points used for computing 
the initial guess point (b) 80x80 grid points used for computing the 
initial guess point (c) 320x320 grid points used for computing the 
initial guess point. 
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Fig. 9. (a) Initial screen of the virtual environment (b) the reaction 
force experienced by over running the particle. 



4. Discussion 

The discreteness of the WE introduced some jagged- 
ness to the haptic sensation. Recall that the method of 
computing the WE used grid points and the smoothness 
of the WE is directly proportional to the grid distance. 
The smaller the grid distance, the smoother is the WE. 
However, the use of force feedback for constraining the 
steering motion of the operator that exceeded the laser 
capability did prove to be a useful method. The reaction 
force from the haptic device made it natural for the 
operator to slow down. In one experiment with the hap- 
tic feedback off, the distance between the operator cur- 
sor and the particle reached up to 3 |Lim. In another 
experiment with the haptic feedback on, the distance 
between the operator cursor and the particle was kept 
under 0.3 |Lim. In order to use this VE in tele-operation, 
much research needs to be done to understand the inter- 
action between the laser beam and the particle and also 
between the particle and the surrounding fluid. 
Validating the model requires deriving a more detailed 
mathematical model for the potential along with exten- 
sive experiments. Also the viscous effect of the sur- 
rounding fluid needs to be taken care of When the par- 
ticle is not a simple sphere, as used in this article, the 
orientation of the particle needs to be considered when 
computing the WE. This means two more variables to 
represent the orientation need to be added to the initial 
three variables (elapsed time, initial velocity in the x 
and J directions). 

When the problem domain is expanded to three 
dimensions, it is presumed that most approaches used 
in two dimensions could be used without major 
changes, with several exceptions. The method for com- 
puting the initial guess point for Eq. (7) cannot be 
extended to three dimensions and needs a new numeri- 
cal method. Also the memory requirement can be sub- 
stantial in three dimensions. The cached data for the 
WE uses 28 MB. An Added dimension (z axis) would 
easily use 1 GB. For a nonspherical geometry one 
would also need to add two more dimensions for the 
orientation. Obviously, a more memory efficient repre- 
sentation of the WE needs to be studied. 



viscous force as a function of a velocity in Eq. (6). This 
would have changed the dynamics described in Eq. (6) 
and Eq. (7). It was not done in this way for simplicity. 
The total processor time to compute this in a batch 
process using an SGI Onyx2 with four processors and 
1 GB of main memory took 71 487.99 s (about 20 h). 



5. Conclusion 

Optical tweezers are a new device that has great 
potential for manipulating nanoscale objects. The non- 
intrusive nature of the device enables nondestructive 
manipulation. To better understand the underlying 
physical nature of the device, a virtual environment that 
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can simulate the physics of the laser beam and particle 
interactions has been designed and implemented. The 
conversion of the system of differential equations to a 
system of difference equations enabled the precise con- 
trol of the transient position of the particle. By putting 
the particle in the maximum force region, the maxi- 
mum movement with the given laser condition was 
achieved. A new concept called WE (Working 
Envelope) and a complex numerical technique to com- 
pute it was also presented. It was shown that WE can 
prevent the user from attempting physically impossible 
movement of the particle. 

The simulation environment could have more value 
when it is used in a tele-operation environment. One 
would have finer position control without ever losing 
the grip. This requires a more sophisticated potential 
model and experiments. 
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